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I. INTRODUCTION 

Mesons are composed of quarks, interacting with each 
other through the exchange of gluons. Under the low 
temperature and density conditions, quarks are confined 
in hadrons and therefore are not observed as free par- 
ticles, which makes their investigation challenging. The 
phenomenon of confinement is closely related to the dy- 
namical chiral symmetry breaking. While the chiral sym- 
metry broken in vacuum is expected to be restored at 
high temperature and density, it is also expected that 
the properties of mesons are affected by this chiral phase 
transition. Thus the investigation of meson properties at 
finite temperature and density is an important part of 
hadron physics. 

Such investigation unfortunately cannot be pursued 
directly in the framework of quantum chromodynamics 
(QCD), the ultimate theory of strong interaction, due 
to its large coupling constant at low energy scale and 
thereby the necessity of dealing non-perturbatively with 
its notoriously complicated structure. The investigations 
are usually performed by using the effective models which 
share the same symmetry properties with QCD such 
as the Nambu-Jona-Lasinio (NJL) model the linear 
sigma model [2] , or by usirig the discretized version of 
QCD, i.e., the lattice QCD Q. 

We employ, in this paper, the NJL model to study the 
nonet meson properties (for reviews on the model, see 
G-g-, 0-01 )• It is not rcnormalizable therefore one has to 
introduce the regularization procedure to handle the di- 
vergences appearing in loop integrals. There are several 
regularization schemes: the three- and four-momentum 
cutoff regularization, the dimensional regularization, the 



Pauli-Villars methods, and so on. We shall use the di- 
mensional regularization (DR) in this paper. In the DR 
method the divergent fermion loop integrals which are 
the functions of the space-time dimensions are regular- 
ized by applying an analytic continuation in the dimen- 
sions variable. Many works are devoted to the studies 
based on it (see e.g., 

We believe that a regularization procedure is an impor- 
tant dynamical part of the NJL model, it gives interaction 
between quarks the size and shapes it, which otherwise is 
point-like in the traditionally used leading order approx- 
imation of the 1/N(. expansion. That is why it is impor- 
tant to chose a proper regularization. Dimensional and 
cutoff regularizations can be considered as corresponding 
to the same size but different shape of the interaction. 
The fact that some physical results depend on the regu- 
larization way implies that the corresponding physics is 
sensitive (probes) not only to the size of the interaction 
but to its shape as well. If so the regularization depen- 
dent properties would be determined by higher energy 
(shorter distances) than those which do not depend on 
the regularization way. 

The purpose of this paper is to study the nonet me- 
son properties such as masses and decay constants in the 
three-flavor NJL model with the DR at finite temperature 
T and chemical potential ^. To this end one has to firstly 
determine the model parameters so that the model de- 
scribes the physical ingredients properly. This has been 
done in the paper [isj . where the present authors have 
found that the various parameter sets reproduce the me- 
son properties nicely at T = and ^ = 0. As our goal is 
to study meson properties at finite T and /i, the paper is 
the straight forward extension of the previous work [l3i | . 

This paper is organized as follows: In Sec.|lTl the three- 
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flavor NJL model and its parameters are introduced. In 
SecHniwe discuss the constituent quark masses and chi- 
ral condensates which are derived through solving the 
gap equations. The results on the meson properties and 
topological susceptibility are shown in Sec. IIVI and |Vl 
Sec. |VT] is devoted to the discussion of the critical tem- 
perature. Summary and discussions are given in Sec. lVIll 
The required calculations of the chiral condensates and 
meson properties are aligned in App.ElandlB] The fitted 
parameters are shown in App. [CJ 



II. NJL MODEL WITH Ua(1) ANOMALY 
A. The model 

The Lagrangian of our three-flavor NJL model is 
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Here the subscripts i,j represent the flavor indices, m 
is the current quark mass matrix diag(m„, m^, TOs), and 
A g are the Gell-Mann matrices in flavor space with Aq = 
•\/2/3 • 1. £4 and Cg are the four- and six-fermion inter- 
actions with the effective coupling constants G and K. 
Cq is introduced to break the Ua(1) symmetry which is 
not realized in the real world. The six-quark vertex in 
the determinant form was firstly discussed by Kobayashi 
and Maskawa and later derived by 't Hooft as an 
instanton-induced quark interaction [15], so this vertex 
is called Kobayashi-Maskawa-'t Hooft term. 

To study the thermal system we evaluate the thermo- 
dynamic potential fl = —InZ/VP, with the partition 
function Z, the volume of the system V and the inverse 
temperature /?(= 1/T). Here we use the imaginary time 
formalism to treat the system at finite T and fj,. Apply- 
ing the mean-field approximation in the DR scheme, we 
obtain the following thermodynamic potential; 
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where 4'i = {ii) represent the chiral condensates and 
Nc{— 3) is the number of colors. Ei — y^p^ + m* ^ is the 
energy of the quasi-particle, and E^^ = Ei zL ij, with the 



In this scheme we regularize the divergent integrals by 
performing the analytic continuation of the space-time 
dimension D to a value less than four, then the integral 
can be written as 
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where Mq is a renormalization scale. One can switch to 
cutoff regularization through the replacement 
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The chiral condensates (jji which are the order param- 
eters of the model are evaluated through solving the gap 
equations, dVt/d4>i = 0. After some algebra, one arrive 
at the following self-consistent equations. 
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with cj„ = (7r//3)(2n + 1), {n = 0, ±1, ±2, • • • ). The tr 
stands for the trace in spinor and color indices. The 
explicit form of the integral is presented in the App. [A) 
To evaluate meson properties, one needs to carry out 
further calculations as discussed in 13| where only the 
case of T, = is presented. We carry out the corre- 
sponding calculations for the case of finite T and fi in 
the App. m 



B. Model parameters 

In the DR, the three-flavor NJL model has seven pa- 
rameters: 

current quark masses rriu, Wd, rus 
four-point coupling G 
six-point coupling K 
dimensions D 
renormalization scale Mq. 

Since the mass difference between to„ and md is small 
as compared to one between rUs and to„ (or rud), we 
take the isospin limit, md = rriu, and test several values, 
niu = 3, 4, 5, 5.5, 6MeV. To fix the other parameters, we 
usually employ the following physical observables [Iq : 



m„ = 138MeV, U = 92MeV, 
niK = 495MeV, m^. = 958MeV. 



(9) 



One more physical quantity is needed to fix all seven 
above listed parameters. Possible candidates are the al- 
ternative meson properties such as the rj meson mass, 
TO,,, topological susceptibility, %, kaon decay constant. 



constituent quark masses m* and a chemical potential fi. /k, etc. We choose the following three cases: 
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Case m 
Case xi70 m^, mK, U, rrirj 
Case xi79 m^, TOK, /tt, 



548MeV 
;^i/4 = i70McV 
X^^^ = 179MeV 



As mentioned in the introduction, the parameters are 
fixed at T, /i = in [l^ ; we show the complete sets of 
them in the App. [C) We do not consider the T and fi 
dependence of these parameters in this paper. In App.O 
we also show the sets in the three-momentum sharp cut- 
off case to make comparison of the two different regu- 
larization schemes. In the cutoff regularization the six 
parameters, mu,md{= 'inu),ms,G, K and cutoff scale A, 
are fixed by using the four physical observables ^ and 
choosing several values of m„ . 



III. CONSTITUENT QUARK MASSES 

Before going through the meson properties, let us dis- 
cuss the results for the constituent quark masses to*. 
This is important because the constituent quark masses 
determine the scale of the system, and they play a cru- 
cial role in the studies of the transition temperature and 
chemical potential. 



m* in the DR scheme 
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FIG. 2: m* in the Case rrir, with various m„. (a) finite T and 
11 = 0. (b) finite p and T = lOMeV. 
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FIG. 1: m* in the Case with m„ = 3MeV. (a) finite T 

and At = 0. (b) finite /i and T = lOMeV. 



In Fig. [U we display the numerical results for the con- 
stituent quark masses which are obtained by solving the 
gap equations Eq.Q as explained in the previous section. 
We note that the absolute values of to* and to* gradually 
decrease with increasing temperature, while at low tem- 
perature they do not depend on the chemical potential 
up to its critical value where they suddenly fall. This is 
a clear signal of the first order phase transition. 

Similar pictures correspond to the different values of 
Tom as it is shown in Fig. [2] where the dependence of 
TO* on T and is given in the Case to^ for fixed values 
of TO„ = 3, 4, 5, 5.5 and 6MeV. These pictures again 
indicate the crossover transition for low /i and the first 
order phase transition for low T. Thus the behavior of the 
first order and crossover transitions holds for the different 
values of to^ in the Case to,,. 

In Fig.[3]there are shown to* for different values of to„ 
as the functions of T or fi in the Case m,^. According 
Eq. ([7]), due to the large contribution of m^, the absolute 
value of TO* is larger than that of to*. We see that the 
qualitative behavior of to* as a function of T (at fixed 
low /i) is similar to that of to*. Yet this is not the case 
with the fj, dependence of to* at fixed low T where there 
appear two gaps (see lower panel). The first gap comes 
from the effect of to* , the second one comes from that of 
TO*. We noted that the critical chemical potential is less 
than the constituent quark mass at T = in the DR for 
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FIG. 3: m* in the Case with various m„. (a) finite T and 
/i = 0. (b) finite and T = lOMeV. 



FIG. 4: m* and m* in the case of cutoff for rriu = 5.5MeV. 
(a) finite T and /i = 0. (b) finite and T = lOMeV. 



2 < £) < 3 at the chiral limit [gj. 

Other three parameter sets, namely Cases Xi70i Xi79 
and Xmi lead qualitatively to similar results. However, it 
is worth mentioning that the Case x?79i although repro- 
duces proper values for the meson properties discussed 
in [l^, shows quite different quantitative behavior. The 
constituent quark masses are one order of magnitude 
larger than the ones in the other parameter sets. This 
leads to unphysical results for the chiral phase transi- 
tion; the transition temperature becomes of the order of 
IGeV which is much larger than the physically expected 
value. This may indicate that the parameters of the Case 
X™ are not suitable to describe the chiral dynamics in 
this model. To explicitly compare these parameter sets, 
we align the numerical results for the constituent quark 
mass, TO* , in Tab. HI There one can clearly see that the 



TABLE I: Constituent quark mass [MeV]. 



TTiu Case Case m,, Case Xi70 Case Xi79 Case xi79 
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-679 


-415 


-420 




6.0 
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values in the Case Xm are considerably larger than the 
ones in the other parameter sets. Since the model predic- 
tions in the Cases Xi70 and X179 are almost the same due 
to the similar values of to* and the model parameters, 
we shall only show the results in the Case Xi70- In what 
follows, therefore, we will focus on the Cases toJ^'^, to^ 
and X170- 



B. m* in the cutoff regularization scheme 

Here we present the results obtained in the widely used 
sharp cutoff scheme to make comparison with our results 
of the DR. 

The numerical results for the constituent quark masses 
TO* and m* for the case to« = 5.5MeV are shown in 
Fig. m It is found that the T and /i dependences are 
similar to the DR case. Fig.[S]shows qualitative difference 
between the cases of to^, = 5.5MeV, and, to„ = 3, 4, 
5MeV. [I dependence of the constituent quark mass to* 
has discontinuity in the case of to„ = 5.5MeV, while to* 
decreases continuously in the cases of to„ = 3, 4, 5MeV. 
Thus, we have the phase transition of the first order in 
the case of mu = 5.5MeV and crossover in the other 
cases. This is a critical qualitative difference between 
the dimensional and cutoff regularizations. 

Fig. [5] displays the corresponding results for to*. We 
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FIG. 5: m* in the case of cutoff, (a) finite T and = 0. (b) 
finite and T = fOMeV. 



FIG. 6: m* in the case of cutoff, (a) finite T and = 0. (b) 
finite and T = fOMeV. 



note that the T dependence of m* resembles each other 
in the two regularizations, which is seen in Figs.|3Ua) and 
mja). However, as seen from Figs. [Sfb) and|6l[b), the ji 
dependence of m* is different. The two gaps in the cutoff 
regularization are smaller and smoother than those in the 
dimensional scheme. 



IV. MESON PROPERTIES 

In this section, we will show the meson properties, 
TTiK, nijj, rrijj/}, in the two regularizations. We 
present the Cases and xno with m„ = 5.5MeV, and 
m^^ with m„ = 3MeV. The reason of these choices is 
that the case with m„ = 5.5MeV is frequently used value 
in the cutoff regularization [1, As the Case m^^ 

does not have solutions at m„ = 5.5MeV, so we select 
the niu = 3MeV case in which the critical temperature 
of the chiral phase transition is close to empirical data 
(see. Sec. EH). 



A. Meson properties in DR 

Here the numerical results of the meson properties are 
presented for the Cases xno, ™r; and mj^^ with some 
discussion on their behavior. 



Fig. [7] shows the results of the meson properties for (a) 
finite T and A* = 0, and (b) finite fi and T = lOMeV 
in the Case xno with m,j = 5.5MeV. As one can ob- 
serve from the upper panel, and stay almost con- 
stant until they cross 2|to* |, then suddenly jump. In the 
similar manner, ttik is nearly constant below the value 
|m* -|-m*|, and it enhances when itik becomes |m* -|-to*|. 
Jtt decreases monotonously and becomes negligibly small 
after m,r exceeds 2|m*|. At low T the dependence of 
meson properties on ^ has sharp discontinuities as seen 
in Fig. [T^b), where the masses TOtt, mx, tti^ suddenly 
get larger and to^/ , /^r falls off. This sharp discontinu- 
ities come from the discontinuous change of the chiral 
condensate, which corresponds to the strong first order 
phase transition. On the other hand in the upper panel 
the meson masses also show jumps as functions of T, al- 
though the chiral condensates change continuously with 
T. This comes from the singular behavior which appears 
particularly in the lower dimensions D < A. Below, we 
will discuss the origin of these discontinuities through the 
analysis of the solution for the pion state. 

In Fig. [5] the function F^{k^) = Re[l - 2A'3n^(fcO)] is 
plotted which determines the pion mass through the con- 
dition J'^(fc°) = (see App. IB ip. Note that all the lines 
diverge at kP — 2|m*|. Due to this divergent behavior 
of J^T^{kP), there is a jump in the dependence of the pion 
mass on T . Namely, we find the unique solution around 
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FIG. 7: Case Xi70 with m„ = 5.5MeV. (a) finite T and ^ = 0. 
(b) finite and T = lOMeV. 
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FIG. 8: T-n (denominator of A^) for three fixed values of T, 
263, 271, 279MeV, in the Case xi70 with m„ = 5.5MeV at 
/X = 0. 
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FIG. 9: Case with m„ = 5.5MeV. (a) finite T and /.t = 0. 
(b) finite /.i and T = lOMeV. 

In Figs. l9l and [TOl vife show the results for the Case 
and mj^^. We observe the similar pictures in these two 
cases; smooth dependence on T at zero ^, and similar 
discontinuities in the dependence on at low T . As is 
seen in Figs.[3I^b) and irUr b). there are two discontinuities. 
The first discontinuity comes from the effect of the gap for 
TO* , and the second discontinuity comes from the one for 
TO*. These gaps are confirmed in Figs. [2jb) and[3jb). It is 
interesting to note that m^i survives at high temperature 
region in Fig. [SJa). This behavior may be attributed to 
the fact that m^/ is smaller than the threshold 2|to*| at 
T = 0. 

We have also calculated meson properties for the cases 
of different m,i, and found that the differences between 
various cases are rather nominal, therefore we showed the 
results for the case of a particular to„. 



fc° ~ 140MeV at T = 263MeV. However, there appear 
three candidates for the solution at T = 271MeV; two of 
them are located near the divergent point and the third is 
larger. We select the larger value, ~ 330MeV, because 
the pion should acquire a larger mass at high T where 
the chiral symmetry is expected to be restored. This is 
the reason why the discontinuous change of m^r at some 
T, is not originated from the sudden change of the chiral 
condensate 0„. In App. IB l1 we give the technical details 



B. Meson properties in cutoff regularization 

Here we present the results for the meson properties in 
the cutoff case. 

Fig. [TT] shows the dependence of the meson properties 
on T and [i. The results in the two regularizations look 
similar to each other. The most noticeable difference is 
in that the curves are rather smooth in the cutoff case. 
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FIG. 10: Case m\f with m„ = 3MeV. (a) finite T and /i = 0. 
(b) finite and T = lOMeV. 
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5.5MeV. (a) finite T and ^ = 0. (b) finite and T = lOMeV. 



There are no discontinuities for low /i in Fig. Illf a) , be- 
cause Eq. (|B11[) does not include divergence for D = A, 
and the solution on the real axis survives at high tem- 
perature region. In Fig. [TlTb) the gap at /i ~ 340MeV 
is smaller than the gap of the DR case. This can be un- 
derstood from the fact that the gap of to* in the cutoff 
scheme is smaller than the one in the DR scheme (see, 



V. TOPOLOGICAL SUSCEPTIBILITY 

The topological susceptibility x is an important quan- 
tity because it is intimately related to the chiral and 
Ua(1) symmetry of QCD. We show the results of x in the 
dimensional and cutoff regularizations in this section. 

Figs. [12] and [12] display the curves for the topological 
susceptibility x as the functions of T or ^ which resemble 
the curves for to* in Figs. [2] and [5] respectively. Thus we 
confirm that the topological susceptibility is influenced 
indeed by the chiral symmetry breaking. 



VI. CRITICAL TEMPERATURE 

We have seen that the chiral condensate and the topo- 
logical susceptibility always show discontinuous changes 
at low temperature in the DR, which indicates the first 



order phase transition. The tendency remains in all cases 
with different parameter sets. On the other hand, the or- 
der of transition is different for the different parameter 
sets in the cutoff case as confirmed in Fig.[5]Jb) and[T3lb). 
This may indicate that the transition pattern crucially 
depends on the regularization procedures. Therefore it is 
interesting to study the values of the critical temperature 
Tc for various cases both in dimensional and cutoff regu- 
larizations. Here we focus on Tc at /i = 0, and make the 
comparison among the two regularization and the lattice 
QCD simulation. 

There are several ways to define the critical temper- 
atures. In this paper, we shall employ two definitions 
given by the maxima of 



dT ' 



and 



dx^ 
dT ■ 



(10) 



The resulting critical temperatures will be denoted by, 
ri"^^ and T^^\ respectively. 

In Tab. HI] the numerical results for the critical tem- 
perature are aligned for various parameter sets. In the 



Case 



with m„ — 4 and 5MeV, we can not find a 



physically meaningful behavior for the phase transition; 
the absolute value of m* becomes larger with increasing 
T. That is why these cases are marked by "none" in the 
table. 

Note that the values of Tc in the DR except for the 
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FIG. 12: Case m^: Topological susceptibility x'^^*- i^) finite 
T and /i = 0. (b) finite and T = lOMeV. 



FIG. 13: Topological susceptibility in the case of cutoff, (a) 
finite T and /X = 0. (b) finite and T = lOMeV. 



TABLE II: Critical temperature 



rriu m,, Xi70 Cutoff 

3.0 184 304 253 146 

4.0 none 336 249 155 

5.0 none 374 243 170 

5.5 - 395 240 184 

6.0 - 419 237 



TABLE III: Critical temperature Ti 



rriu itlj) Xno Cutoff 

3.0 188 309 259 152 

4.0 none 341 255 161 

5.0 none 379 250 175 

5.5 - 401 246 189 

6.0 - 425 244 



Case mj^ are larger than ones in the cutoff case. The 
Case m}^ at m„ = 3MeV leads almost to the same value 
of Tc as in the cutoff case at m„ = 5.5MeV which is 
comparable with the result in the lattice QCD simula- 
tion, 150 - 200MeV We see that becomes larger 
with increasing m„ in the Case while does not 
change drastically with m„ in the Case Xito- This can 
be explained by the values of to*. 



VII. SUMMARY AND DISCUSSIONS 

In this paper we studied the nonet meson characteris- 
tics, m^, tok, TOr,, TOr,', and x finite temperature 
and chemical potential in the NJL model with DR. Many 
of meson properties show reasonable behavior, which is 



similar to one obtained in the model with the cutoff rcg- 
ularization. 

We examined the behavior of the constituent quark 
masses, to* and to*, in the two regularizations. T de- 
pendence of the constituent quark masses at low /i does 
not depend essentially on the way of regularization. How- 
ever, /X dependence of the constituent quark masses at low 
T depends on the regularization way around the critical 
chemical potential. In the DR, the phase transition at 
low T is always of the first order for various choices of rriu- 
This tendency is consistent with the current consensus 
about the transition to be of the first order in the chiral 
limit [l^. On the other hand in the cutoff case, the gap 
at the critical chemical potential becomes smaller with 
decreasing value of to„, and it eventually vanishes result- 
ing in crossover transition (see, Fig. [5]). This happens 
due to the fact that the larger cutoff A diminishes the 
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coupling strength which leads to the smoother change of 
the order parameters. Decreasing the coupling strength 
weakens nonperturbative effects. This explains the effect 
of the UV-cutoff which is also seen in the work 120^] as 
the effect of the cutoff in the temporal direction. 
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Using the obtained values of the constituent quark 
masses, we evaluated meson properties at finite T and/or 
/i in the two regularizations. T dependence of the meson 
masses at low ji are similar to each other except around 
mp ~ I TO* + TO*|. Since the singularity of the pion self- 
energy np(fc^) appears at fc^ = (m* + ■m*Y for £> < 4, 
the behavior of meson masses becomes discontinuous in 
the DR. n dependence of the meson masses at low T also 
does not essentially depend on the way of regularization 
except at the region of the critical chemical potential. 
Since the gap of m* in the DR case is larger than one in 
the cutoff case, the gap of the meson mass becomes large. 
The behavior of /^r as the functions of T or in the DR 
case resembles one in the cutoff case. We also evaluated 
the topological susceptibility x as the functions of T or 
(jL in the two regularizations. The behavior of x is similar 
to one of TO* . 



In the DR, we examined the three parameter set cases: 
Xi70, m^i and . Parameters of each set are fixed by 
fitting physical quantities at T, /i = [l^ . Most of the 
properties concerning meson masses and decay constants 
are similar to each other in the three cases. However, the 
large difference is seen in the values of the critical tem- 
perature and of the critical chemical potential. One of 
the sources of this difference may be the strength of cou- 
plings, i.e., increasing the values of the coupling strength 
decreases the values of the critical temperature and the 
critical chemical potential. We actually calculated the 
critical temperatures, Tc'^^ and Tc^\ In the Case rn}^ 
at TO„ — 3MeV, the critical temperature is almost the 
same as in the cutoff case at m„ — 5.5MeV which is 
comparable to the result in the lattice QCD simulation. 



The lessons we learned in this paper are twofold: (i) 
similarities in the predicted meson properties and (ii) re- 
markable differences seen in the order of the chiral phase 
transition at finite chemical potential. Regarding the 
point (i), we have been convinced that the NJL model 
predictions do not depend drastically on whether one uses 
the dimensional or cutoff regularization procedure. This 
may be the regularization independent aspect of the NJL 
model, which is intriguing problem to be studied in more 
detail in future. As to the point (ii), the results of the 
paper suggest that the critical point in the phase diagram 
on T — /i plane shows different behavior. We think that 
it is interesting to study the behavior of the chiral phase 
transition and the phase diagrams of the model in more 
detail, especially in context of how the location of the 
critical point depends on the parameters in the different 
regularization schemes. 



Appendix A: Chiral condensate 

By taking the trace with respect to the Dirac spinor 
indices in Eq. (jSl), we get the chiral condensate 



= trl • to" 



T 



(Al) 



where trf = ■ 2^/"^. After performing the frequency 
summation, we arrive at the expression 



trl • to' 



jD-1 



P -1 



(27r)^-i 2E, 



, (A2) 



with the Fermi-Dirac distribution f{E) given by 



1 



1 + e^^. 



(A3) 



Appendix B: Meson properties 

Here we evaluate meson masses, decay constants and 
topological susceptibility. More detailed, self-contained 
calculations are presented in the review paper 



1. TT and K masses 

Meson masses are obtained through examining the pole 
structure of their propagators. Using the random-phase 
approximation and the 1 /A'^^expansion, we get the meson 
propagators of the form d, Q 



Ap(A:2) 



2K^ 



1 - 2X„Hp(fc2; 



+ 0{N-^), (Bl) 



where the index P(= vr, K) denotes the meson species 
and a labels the isospin channel. The pole position is 
determined by the equation 



1 - 2KJ1y>{x^) = 



(B2) 



where is regarded as the square of a meson mass, i.e.. 



TOi or ml 



The effective couplings, K^, of mesons become 
ifg = G - iif 0,, for ttO, 
K^ = G-^K4>u, for K°,K°. 



(B3) 
(B4) 
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The self-energy, Hp, of the propagator is given by 

= / l^^D^''\ri^^^S\p + kl2)^,TlS^{p - A/2)] , 

(B5) 

where the trace runs over flavor, spinor and color indices. 
The flavor S'C/(3) matrices are: T3 = A3 for 7r°, Te = 
(A6 + iA7)/V2 for and Tj for K". Then the self-energy 
for 7r°, K° and K'' can be written as 

n,(fc2 = ml) = 2nr(fc' = ml), (B6) 
nK(fc' = ml) = 2nr(fc' = ml), (B7) 

where ng"'(fc^) is defined by 
n?(*:")=y'-|^tr[75S'(p + fc/2h,,Si(p-*V2)] 



with 



trl 



(B9) 



The integral above becomes 



- /z(2^)^^ 1 1 ^. 2K. dI ^^^'\ 



^0 (.i^j 



D 



°-[f{E+)-f{En-f{El) + f{Er)] 



(BIO) 



where Of^ = {Ei ± Ej)'^ ~ k'^ and = m* ^ - m* ^ - fc^. 
In practice we first solve the gap equations to get the 
chiral condensates 0; and constituent quark masses m*, 
then we scan the solution of Eq. (jB2l) in evaluating TOtt 
and mx- In the actual numerical calculations, we choose 
the rest frame A:^ = (A;°, 0). 

It should be noted that Itj as a function of [k = 0) 
has a singularity. Since the pion self-energy 11^ consists of 
the function it is divergent at the singularity, n^(fc^) 
has a cut in the region fc" > 2|to*| on the real axis. One 
should move the other Riemann sheet to calculate IiTr{k^) 

(v) . 

at the cut. The vacuum part /„„ is written as 



N, 



where 



(27r)^/: 



Luu {k 



,r(2. 



) dxL^l^-\k% (Bll) 



(B12) 



k^xil 



The function /„ti(^^) is divergent at fc^ = 4m* ^ for D < 
4. Then Il^{k^) has the singularity at = 2|rn*| which 



is confirmed in Fig. [8] where D = 2.47 (Case xno with 
m„ = 5.5MeV). The divergence appears at the resonance 
position. It does not mean the existence of a bound state 
except for k^ = 2|m* |. In Fig. [5] three solutions are 
observed for J"^ = at T = 271 and 279MeV. Two of the 
solutions are located near the singularity at fc° = 2|r7i*|. 
We take the solution with the larger to express the 
pion mass. Thus the meson masses jump as functions of 
T near the critical temperature. 

As is shown in Fig. HI the larger solution disappears 
for T > 279MeV. In this case we determine the meson 
mass as fc° at the maximum of J-'^. On the contrary, the 
reliable solution always exists and the smaller two solu- 
tions are degenerate above the critical chemical potential 
at T = lOMeV. 



2. TT and K decay constants 

Here, we display the equations determining /p, the 
pion and kaon decay constants, with P(= tt, K). 



ikfj.fpdal3 



{2tt)D 



tr 



T 

7/<75 ^ S'gpgg (0)75TjS'-' 



(B13) 



where the renormalization scale Mq is introduced so 
that the decay constant has a correct mass dimension, 
dim(/p) = 1. The trace is with respect to the flavor, 
spinor and color indices. The meson-to-quark-quark cou- 
pling, gpqq, is defined by 



4_^9IIp(fc2) 

9/c2 ■ 



Eg. dBlSp can be rewritten in the following form 

/2 = -M^-''gpqqll/k\ 

where /^r = TO*/,ju(/c2)fc^ and 

d^P (m* - ml)p^ - mlkf, 
i{2iT)^ (p'^ - ml^)[{p - k)"^ ~ ml"^ 



(B14) 



(B15) 



trl 



By inserting the values of the constituent quark masses 
m* and the renormalization scale Mq, one can easily eval- 
uate the decay constants. 



3. rj and rj' mesons 

The equation which determines 77 meson masses be- 
comes 2x2 matrix due to the well known Ua(1) anomaly. 
In the leading order of the 1 /N^ expansion, the propaga- 
tor of the rj — rj' system becomes 



A+(/fc2) = 2K+[1 - 2K+n(fc2)]-i, 



(B16) 
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where K"*" and 11 are given by 



with 



n 



K[)n = G 



and 













Hob 




HsB 



i^BB = 



6 



K. 



%/2 



BO 



K{4>s - <t>u), 



noo(fc2)--[2nr(fc2) + nr(fc2)], 



n88(fc2)^_[n^"(fc2) + 2n^^(fc2 



no8(fc') -n8o(fc2) 



2^2 



(B17) 
(B18) 



TABLE IV: Case m^°. 





TTls 


G 


K 






B 


3.0 


84.9 


-0.0195 


9.02 X 10" 


7 


118 


2.29 


4.0 


118 


-0.0174 


9.17 X 10" 


7 


113 


2.38 


5.0 


156 


-0.0162 


9.49 X 10" 


7 


108 


2.47 






TARTjE V 










niu 


nis 


G 


K 




Mo 


D 


3.0 


79.0 


-0.0130 


2.29 X 10" 


7 


107 


2.37 


4.0 


106 


-0.00748 


8.26 X 10" 


8 


92.0 


2.52 


5.0 


134 


-0.00357 


1.99 X 10" 


8 


73.2 


2.69 


5.5 


147 


-0.00231 


8.40 X 10" 


9 


62.4 


2.77 


6.0 


162 


-0.00142 


3.23 X 10" 


9 


50.9 


2.87 


TABLE VI: Case Xito- 


niu 


nis 


G 


K 




Mo 


D 


3.0 


77.1 


-0.0168 


2.23 X 10" 


-7 


120 


2.28 


4.0 


106 


-0.0143 


2.11 X 10" 


-7 


116 


2.36 


5.0 


134 


-0.0119 


1.80 X 10" 


-7 


112 


2.43 


5.5 


150 


-0.0109 


1.62 X 10" 


-7 


110 


2.47 


6.0 


166 


-0.00992 


1.48 X 10" 


-7 


109 


2.50 



The condition which determines the 77 and 77' masses is, 

det[l - 2K+n{x^)] = 0. (B19) 

Using the constituent quark masses derived from the gap 
equations, we search for the appropriate solutions for 
and m,i' . 

4. Topological susceptibility 

The topological susceptibility x is given by the corre- 
lation function between the topological charge densities 
Q{x) at different points. 



connected 



(B20) 



here the charge density is defined as 



Qix) = = 2XIm[detg"(l - 75)9], (B21) 

with g, the strong coupling constant of QCD, and i^"^, 
the gluon field strength. The explicit form of x is given 



below m 



X = 



M, 



D-i ^" 



20s 
ml 



V6 



(20s + <^„)(noo(o),no8(o)) 



V3 



(0, (no8(o),n88(o)) I A+(o) 



V J ^ (0^ noo(o) 



no8(o) 
nB8(0) 



(B22) 



After solving the gap equations, one can numerically cal- 
culate this quantity by inserting the values of (pi and m* . 



Appendix C: Parameter sets 

Here we present the parameter sets obtained in [Tl]. 
In the Case tti^ and Xi79j there exist two sets of param- 
eters, and we distinguish between these cases by using 
the superscript LD (lower dimension) and HD (higher 
dimension) . 

In Tab. IIXI we also present the parameter sets of the 
Case Cutoff. In the cutoff regularization, the number of 
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TABLE VII: Case xi79 TABLE IX: Case Cutoff. 



ruu 


rus 


G 


K 


Mo 


D 




rUs 


GA^ 




A 


3.0 


78.0 


-0.0170 


2.70 X lO"'' 


120 


2.28 


3.0 


89.5 


1.55 


8.34 


960 


4.0 


106 


-0.0144 


2.52 X lO"'' 


115 


2.36 


4.0 


110 


1.60 


8.38 


797 


5.0 


136 


-0.0120 


2.11 X 10"'^ 


111 


2.44 


5.0 


128 


1.71 


8.77 


682 


5.5 


152 


-0.0109 


1.90 X lO"'^ 


110 


2.47 


5.5 


136 


1.81 


9.17 


630 


6.0 


168 


-0.00995 


1.70 X lO"'' 


108 


2.50 


5.87 


139 


2.09 


10.1 


580 



TABLE VIII: Case xiii 
TUu rUs G K Mo D 

3.0 74.8 —4.24 x 10~^ 1.38 x 10~^^ 2.91 3.22 parameters is five, and these are fixed by using the four 
4.0 100 -6.91 X 10"^ 3.60 x 10"^^ 1.15 x 10"^^ 3.90 physical observables m^r, /tt, mx, m,,' for several values 

— ofm„. 
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